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Preface 


With  the  current  emphasis  of  fighter  development  on 
sustained  high  performance,  the  chances  for  aircraft  losses 
due  to  the  pilot  entering  a flight  condition  from  which  the 
aircraft  will  not  recover  become  very  great.  Ideally,  design 
of  the  aircraft  includes  aerodynamic  characteristics  or 
augmentation  which  make  the  aircraft  free  from  the  poss- 
bility  of  loss  of  control.  For  the  newer  generation  of 
fighters  which  have  a high  degree  of  augmentation  already 
actively  part  of  the  flight  control  system,  the  possibility 
of  a departure  controller  being  added  offers  a relatively 
inexpensive  means  of  improving  both  aircraft  performance 
and  flight  safety.  It  was  my  intention  to  determine  for  the 
YF-16  where  in  the  alpha-beta  control  boundary  loss  of  con- 
trol might  be  a serious  problem.  I then  hoped  to  be  able 
to  put  the  uncontrolled  motions  resulting  from  pitch  and 
yaw  perturbations  into  small  enough  categories  so  that  a 
simple  controller  could  prevent  the  motions  by  monitoring 
a minimum  number  of  variables. 

I wish  to  thank  my  thesis  advisor,  Dr.  Robert  A. 

Calico,  of  the  Department  of  Aeronautics  and  Astronautics 
of  the  Air  Force  Institute  of  Technology.  He  was  instru- 
mental in  helping  me  tackle  this  project  with  a degree  of 
organization  and  insight  that  I could  not  have  attained 
without  guidance. 

I also  would  like  to  thank  Tom  Chord  and  the  other 
members  of  the  Aircraft  Group,  Control  Dynamics  Branch, 
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Flight  Control  Division  of  the  Air  Force  Flight  Dynamics 
Labratory  who  helped  me  become  familiar  with  the  computer 
simulation  program  which  formed  the  basis  of  this  study. 

Host  importantly,  I wish  to  thank  my  wife,  Marilyn, 
for  her  support  and  devotion  during  this  very  hectic  time 
in  our  lives.  This  thesis  is  dedicated  to  her  and  to  our 
new  son,  Ryan. 
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Abstract 

A, 

A study  was  made  of  the  theoretical  departure  modes  of 
the  YF-16  due  t o pitch  and  vaw  perturbations  from  asymmetri- 
cal rectilinear  flight.  An  alpha-beta  control  boundary  was 
developed,  ami  within  that  boundary  perturbations  of  .’0  dee,/ 
sec  and  30  deg/ see  were  introduced.  The  areas  of  uncontrolled 
motions  were  mapped  out  and  the  motions  were  categorized. 

Three  types  of  motions  were  identified.-  erect  spins,  inver- 
ted spins,  and  rolling  departures.  For  yaw  perturbations  onlv 
a simple  controller  was  modelled  which  had  as  its  inputs 
alpha,  beta,  and  yaw  rate.  The  controller  prevent  oil  the  air- 
craft from  departing  controlled  flight.  Time  traces  ot  the 
various  departure  modes  and  the  effects  ot  including  a depar- 
ture controller  are  presented  in  the  appendices. 


INVESTIGATION  OF  THE  YF-16  IN  HIGH  ANGLE 
OF  ATTACK  ASYMMETRIC  FLIGHT 

I . Introduction 

Background 

With  the  development  of  a new  generation  of  fighter 
aircraft  capable  of  ever  increasing  performance,  there  is  a 
renewed  interest  in  efforts  to  eliminate  susceptibility  to 
loss  of  control  or  spin  entry.  Attempts  have  been  made 
to  provide  design  guidelines  so  that  new  aircraft  will  be 
far  more  departure-resistant.  These  guidelines  have  their 
limitations,  however.  Invariably  a new  fighter,  once  opera- 
tional, will  have  some  portion  of  its  flight  envelope  where 
through  flight  tests  or  analytical  studies,  spin  susceptibil- 
ity is  either  experienced  or  predicted. 

The  new  fighters,  because  of  the  high  degree  of  flight 
control  system  augmentation,  have  the  possibility  of  rela- 
tively inexpensive  modifications  being  done  to  correct  or 
enhance  overall  performance.  The  task  must  be,  then,  to 
identify  potentially  dangerous  flight  regimes  and  to  analyze 
how  a flight  control  system  (FCS)  can  be  used  to  help. 

Purpose 

The  purpose  of  this  study  was  to  investigate  the  theo- 
retical departure  potential  of  the  YF-16  due  to  the  influ- 
ence of  pitch  and  yaw  perturbations.  The  resulting  uncon- 
trolled motion  was  to  be  categorized  as  much  as  possible. 


Trimmed  asymmetric  flight  was  used  as  the  initial  condition. 
The  portions  of  the  a - 6 envelope  where  departures  might  occur 
were  to  be  identified.  If  results  showed  regions  where  a 
small  set  of  variables  were  departure  producing,  then  a de- 
parture controller  which  was  a function  of  those  variables 
could  be  modelled  as  part  of  the  FCS . 

Model 

The  YF-16  was  selected  as  the  model  for  this  study 
because  it  represents  the  newest  concept  in  fighter  design, 
that  of  static  instability.  This  feature  makes  the  aircraft 
both  highly  maneuverable  and  highly  augmented.  During  flight 
tests,  the  YF-16  had  certain  flight  conditions  identified  as 
potentially  hazardous.  Two  departures  and  one  actual  spin 
were  encountered.  The  aerodynamic  data  available  was  also 
of  interest,  as  General  Dynamics  had  done  studies  to  incor- 
porate the  data  obtained  from  the  spin  and  the  departures 
into  a new  data  package  with  simplified  functional  depend- 
encies from  the  previous  data  package.  The  physical  char- 
acteristics of  the  YF-16  that  were  used  for  this  study  are 
listed  in  Appendix  A. 

Data  Package 

The  aerodynamic  data  for  the  YF-16  was  obtained  from 
a study  done  by  General  Dynamics  in  April  1977.  This  data 
was  used  because  extensive  work  had  been  done  to  update  all 
aerodynamic  coefficients  based  on  the  results  of  the  flight 
test  evaluation  at  Edwards  AFB  in  1975.  This  revised  data 
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package  not  only  modelled  actual  aircraft  flight  tests  more 
closely,  but  was  expanded  to  provide  data  up  to  '120  degrees 
angle  of  attack  and  up  to  degrees  sideslip  angle. 

A significant  advantage  of  the  expanded  data  package 
is  that  the  functional  dependencies  of  the  aerodynamic  coef- 
ficients are  reduced,  thereby  simplifying  the  calculations 
necessary  for  digital  computer  simulation.  Also,  the  equa- 
tions developed  for  the  linear  theory  portion  of  this  study 
are  more  compact  than  would  have  been  possible  with  the  ex- 
isting YF-16  data  package. 

The  functional  dependencies  of  the  aerodynamic  coef- 
ficients are  as  follows: 
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The  coefficients  obtained  from  General  Dynamics  and  the  inde- 
pendent variables  for  the  coefficients  are  listed  below: 


(1)  c*  - cw(“-  5s>  + V"’  H <7) 

(2)  C = C (a)  + C (a)  6 + C («)  & (8) 

y y*  y^  a y&  r 

u r 

<3>  Cz  - C <„.  <„>  + C2  („)  « (9) 

basic  q 2V 

I 

3 

l><ll*<M^^^*******l**,*l**lll*||*M— — I -ad 
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(4)  Cx  = Cj_  (a,  3)  + Cx  (a)  6fl  + Cj_  (a)  6r 

basic  <5  6 

a r 

+ C (a)  pb  + C (a)  rb 

Lp  T7  lr  2V 


<5>  C“  " S.asic^'^  + S<0)  i 


(10) 

(11) 


(6)  Cn  = C (a,  0)  + Cn  (a)  + Cn  (a)  6r 

Dasic  6a  6r 


+ C (a)  pb  + C (a)  rb 
np  2V  nr  5V 


(12) 


The  effects  of  all  other  control  surfaces,  such  as  the  auto- 
matically scheduled  leading  edge  flaps  and  the  differential 
horizontal  tail,  were  incorporated  into  the  aerodynamic 
coefficients  for  the  three  primary  control  surfaces  by 

General  Dynamics.  The  range  of  values  for  a,  0,  and  6 are: 

s 

(1)  a:  from  -120  to  +120°  in  5 increments;  49  values 

(2)  0:  from  -40  to  + 40°  in  5 increments;  9 values 

(3)  6 : -25,  -10,  0,  +10,  +25;  5 values 

s 

In  addition  to  these  coefficients,  the  following  coef- 
ficients were  derived  using  finite  differences  techniques 
(Ref  6): 


(1)  C (a,  6S),  C (c,  6S),  Cm  (c.  « ) 

a a a 

(2)  C,  (a,  0),  C,  (a,  0),  C (a,  0),  C (a,  0) 

Act  A0  na  n0 

(3)  C (a.  Ss),  C (a,  y,  C (a,  4,) 


These  coefficients  were  necessary  for  the  perturbation  equa- 
tions used  in  the  eigenvalue  problem.  All  of  the  data  ob- 
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tained  from  General  Dynamics  was  for  stability  axes.  The 
data  was  converted  to  body  axes  since  the  non-linear  computer 
program  was  written  for  body  axes. 

Method 

From  the  flight  tests  the  areas  of  potential  departure 
problems  in  terms  of  a and  0 were  fairly  well  known. 

However,  the  first  step  in  the  analysis  was  to  define  an 
a-0  control  boundary.  This  was  accomplished  by  use  of  a 
computer  program  developed  by  Hawkins  (Ref  4)  which  trimmed 
the  aircraft  for  level,  unaccelerated  flight.  The  values  of 

o o 

a ranged  from  0 to  29,  which  was  the  maximum  attainable 
angle  of  attack  for  level  flight  for  the  YF-16  during  its 
flight  tests.  The  values  of  6 ranged  from  0 to  40.  The 
aileron-rudder  interconnect  (ARI) , which  serves  to  help 
coordinate  turns,  was  not  modelled  for  the  trim  program. 

This  resulted  in  slightly  more  control  authority  for  cross 
controlling  the  aircraft  than  was  available  for  the  actual 
aircraft  during  flight  tests.  The  solution  to  the  trim  equa- 
tions required  that  density  be  constant.  Therefore,  the 
a-0  boundary  was  arbitrarily  developed  for  sea  level. 

The  second  step  was  to  test  each  point  in  the  control 
limits  boundary  for  stability.  The  linearized  equations  of 
motion  were  used  in  this  analysis.  It  was  necessary  to  include 
the  longitudinal  flight  control  system  when  testing  for 
stability  since  the  aircraft  was  statically  unstable  for 
practically  all  of  the  a-P  region  that  was  considered.  A 
state  space  approach  was  used  to  obtain  the  eigenvalues  and 
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eigenvectors  of  the  system.  The  stability  of  all  trim  points 
from  the  trim  program  was  checked. 

The  third  step  was  to  test  the  behavior  of  the  aircraft 
trimmed  for  level  flight  and  to  determine  what,  if  any, 
adverse  flight  conditions  might  be  encountered  when  the  air- 
craft was  perturbed  from  an  equilibrium  trim  position.  Since 
this  study  was  interested  in  determining  what  severe  situa- 
tions might  be  encountered,  the  sizes  of  the  perturbations 
were  increased  until  the  aircraft  experienced  such  conditions. 
The  equilibrium  flight  conditions  from  which  the  aircraft  was 
disturbed  represented  slow  speed,  level,  or  near  level,  unac- 
celerated flight  with  landing  gear  and  trailing  edge  flaps 
up . 

Finally,  for  the  uncontrolled  flight  conditions  that 
were  encountered,  an  analysis  was  made  to  determine  if  these 
conditions  could  be  grouped  into  any  general  categories.  It 
was  desired  to  determine  a minimum  number  of  independent 
variables  that  could  describe  regions  where  loss  of  control 
was  imminent,  and  to  use  these  variables  to  model  a depar- 
ture controller. 
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II.  Linear  Analysis 


Development  of  the  Linear  Equations  of  Motion 

The  system  of  equations  used  for  both  the  trim  analysis 
and  the  stability  analysis  is  that  found  in  Reference  3:149. 
These  equations  assume  the  earth  as  an  inertial  reference 
frame,  an  atmosphere  at  rest  with  respect  to  the  earth,  and 
a rigid  aircraft.  The  equations  in  body  axes  notation  are: 


X 

- mg 

sinO 

= 

m(u 

+ 

qw  - rv) 

(13) 

Y + mg 

cosO 

sin4> 

= 

m(v 

+ 

ru  - pw) 

(14) 

Z + mg 

COS0 

costf> 

= 

m(w 

+ 

pv  - qu) 

(15) 

L 

* 

XxP 

- 

Zzx<*  + pq)  - (Iy  ' 

!z)qr 

(16) 

M 

xr 

V 

- 

Wr2-  P2>  - (Iz  - 

Ix)rp 

(17) 

N 

*= 

V 

- 

WP  * qr)  ' (Ix  “ 

!y)pq 

(18) 

♦ 

= 

P + 

q 

sin^  tanO  + r cos^ 

tanO 

(19) 

0 * q cos4>  - r sin4>  (20) 

^ = (q  sinij>  + r cos4')/cos0  (21) 

This  is  the  system  of  three  force  equations,  three  moment 
equations,  and  three  kinematic  equations  which  had  to  be 
linearized  and  equilibrium  points  determined. 

Trim  Equations 

It  was  desired  to  have  the  aircraft  trimmed  for  asym- 
metrical rectilinear  flight.  Thus,  ae,  Be>  0e,  4>e , and  m 
were  constants.  Equilibrium  was  defined  such  that  p,  q,  r, 
u,  v,  w,  p,  q,  and  r were  zero.  Also,  the  aircraft  was 
assumed  to  be  symmetric  about  the  xz  plane,  and  thrust  was 
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assumed  to  act  along  the  x-axis.  For  these  assumptions,  equa- 
tions (13)  through  (18)  reduce  respectively  to  the  following: 


Tp  + Xp  - mg  sinOe  = 


Yg  + mg  cos 6 e sin^ 

Z + mg  cos 
e ° e 


sint)e 

= 0 

(22) 

sin*e 

= 0 

(23) 

cos 

e 

= 0 

(24) 

Le 

= 0 

(25) 

Me 

= 0 

(26) 

Ne 

= 0 

(27) 

Non-dimensionalizing  equations  (22)  through  (27)  according  to 
Reference  3 yields: 


CT  + Cx  ' Cw  sin0e 
e e e 


C + Cw  cosOe  sin^e 
* e e 


C + C cosO  cosJi 
e ' e 

e e 


m 


n 


0 

0 

0 

0 

0 

0 


(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


The  force  and  moment  aerodynamic  coefficients  are  expressed 
in  the  following  functional  forms  for  a,  8,  and  being  con- 
stants : 


Cx<“e’  6s  > 


Cy  = C-  (ae)6a  + C (ae)6r 


(34) 

(35) 


8 


(36) 


C (a  .6  ) 

z e se 


" cl(^e'?e)  + C1  (ae)  6a  + Cl„  (ae>  6r 


m 


'n 


e 


Cn<V*e>  + Cn » <"e>  6a  + Cn<  <*e>  5r 
6a  5r 


(37) 

(38) 

(39) 


Substituting  equations  (34)  through  (39)  into  equations  (28) 
through  (33)  gives: 


z w 


Ja 


n 


n, 


C sinO 

we  e 

= 0 

(40) 

cosOe  sin^e 

= 0 

(41) 

COS  0 COS* 

C C 

= 0 

(42) 

6 + C,  5 

a Iff 

6r 

= 0 

(43) 

C 

m 

= 0 

(44) 

6 + C 5 

a n r r 

= 0 

(45) 

The  unknowns  in  the  above  equations  are  6 , 6 , C , and  C_  . 

M a r wg  Te 

Since  Sg  can  not  be  explicitly  expressed  in  terms  of  the 
aerodynamic  coefficients  due  to  it  being  one  of  the  inde- 
pendent variables  for  the  coefficients,  is  obtained  by  an 
iterative  technique.  Equation  (44)  is  satisfied  by  scanning 
the  tabulated  values  for  Cm,  which  is  listed  for  a given  a 
and  6 . The  value  of  6g  is  found  where  for  a given  a,  Cm  is 
zero.  Once  5g  is  known,  the  other  unknowns  can  be  found 
from  equations  (40)  through  (43)  and  equation  (45) . Solving 
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equations  (43)  and  (45)  simultaneously  gives 


ClCn  Cn 

6a  dr 


C,  C - C,  C 

l5  n6a  14  nvS 

r a a r 


and 


(46) 


4.  ' ' cn  , 4r  ’ cn 


(47) 


n 


5 


Dividing  equations  (41)  by  equation  (42)  gives: 


<f>  * tan" C <S  + C 6 , 

e / Vs  a y.  r' 

lSa  '•V 


(48) 


Equation  (42)  can  now  be  solved  for  C 


w 


w 


- C. 


COSO  COS^e 


(49) 


Equation  (40)  gives 


CT  = Cw  sinOe  - Cx 
e e 


(50) 


Also,  V and  T are  calculated  from  the  following  relations: 
e e ° 


2W 

^7^ 


(51) 


*•>.  Ve  S C1 


(52) 


A computer  program  written  by  Hawkins  and  revised  for  the 
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differences  in  functional  dependency  of  the  data  in  this 
study  was  used  to  calculate  all  of  the  above  unknowns.  The 
iterative  approach  for  finding  As  was  a modification  for  this 
study.  The  data  from  the  trim  program  was  then  prepared  for 
use  in  the  eigenvalue  problem.  In  addition  to  this  trim  pro- 
gram, a comparison  of  values  obtained  from  the  non-linear 
program  (Appendix  for  the  same  equilibrium  conditions  was 
made.  Agreement  between  the  two  approaches  was  good  with  the 
exception  that  due  to  the  ARI  being  in  the  non-linear  program, 
slight  variations  of  3 and  $ were  noted. 

Deve lopment  of  the  Unaugmented  Tertvirbat  ion  Kquat  ions 

The  approach  to  the  stability  analysis  was  to  use  the 

linearized  equations  of  motion.  As  in  Hawkins'  study  (Ref  4), 

small  perturbation  theory  was  used  to  linearize  the  equations. 

An  equilibrium  of  asymmetric  rectilinear  flight  was  used. 

This  requires  that  p = q = r =0.  It  is  assumed  that  p 

e e e e 

is  constant,  although  slight  climbs  or  dives  are  permitted. 

The  linearization  of  equations  (13)  through  (31)  using  small 
perturbation  theory  results  in  the  following: 


AX  - mg 

cosO 

AO  = m(Au  + w Aq  - v Ar) 

(53) 

e 

e e 

AY  + mg 

cosO 

e 

cos*1  A4>  - mg  sinO  s i n t AO 
e e e 

* m(Av  + u^  Ar  - wfc,  Ap) 

(54) 

AZ  - mg 

cost1  „ 
e 

sin^  A|  - mg  s inO 0 cospo  AO 

« m(Aw  + v Ap  - u Aq) 
e e 

(55) 

AL  * Ix  Ap  - I2X  Ar 

(56) 

AM  - Iy  Aq  (57) 
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IP 


AN  = Iz 

Ar  - Izx  Ap 

(58) 

A 4>  = Ap 

+ Aq  sin-,*1  tant3  + Ar  cos*  tanO 
. e e e e 

(59) 

A0  = Aq 

cost|)g  - Ar  sin^g 

(60) 

Eq . 13  is  not  used  since  it  is  uncoupled  from  the  other  equa- 
tions. It  was  desired  to  use  AV,  Aa,  and  AS  as  state  vari- 
ables for  the  eigenvalue  problem  rather  than  Au,  Av,  and  Aw. 


Therefore , 

the  following  relationships 

are  used  (Ref  3:117): 

u = V 

cosa  cosS 

(61) 

> 

ii 

> 

sinS 

(62) 

> 

ii 

> 

sina  cosS 

(63) 

Lineariza tion  of  eqs.  (61)  through  (63)  results,  respectively, 
in: 


Au  = cosae  cosSg  AV  - Vg  sinag  cosBe  A a 

- Vg  cosag  sinBg  AS  (64) 


Av  = 

sinSg 

AV  + Vg  cosSg 

AS 

(65) 

Aw  = 

sinae 

cosSg  AV  + Ve 

cosa.  cosSQ  Aa 

c t: 

- Vg  sincig  sinSg  AS  (66) 


The  aerodynamic  forces  and  moments  were  assumed  to  have  the 
following  functional  dependencies : 


AX  = AX(V,a,q, 5g) 
AY  = AY(B,Sa,5r) 
AZ  = AZ(V, a, q , 5g) 
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(67) 

(68) 
(69) 


r 


AL  = AL(a, B.p.r, 6a, 6r) 
AM  = AM  ( a , q , 6 g ) 

AN  = AN(a, P,p,r, 6a, 6r) 


(70) 

(71) 

(72) 


Linear  expansions  of  eqs . (67)  through  (72)  results  in: 


AX  = Xv  AV  + Xa  Aa  + X Aq  + Xg  A6( 


AY  = Ya  Ap  + Yx  A 6 + Y , A<S 

6 6a  a 6r  r 

AZ  = Zv  AV  + Za  Aa  + Z Aq  + A6, 


(73) 

(74) 

(75) 


AL  = La  Aa  + L„  AP  + L Ap  + Lr  Ar  4-  Lg  A6  + L6  AS  (76) 

P a a r 


AM  = M Aa  + M Aq  + Mc  AS 
a q ^ 5 s 

s 


(77) 


AN  = Na  Aa  + Ng  Ap  + N Ap  + Nr  Ar  + Nfi  A6&  + Ng  A6r(78) 

p a r 

Substituting  eqs.  (64)  through  (66)  and  similar  expressions 
for  Au,  Av,  and  Aw  into  eqs.  (53)  through  (60)  gives: 


Xy  AV  + Xa  Aa  + X_  Aq  + X^  A5g  - mg  cos0e  A0 

s 

= m(cosae  cosPe  AV  - Vfi  sinae  cosPe  Aa 


- V„  cosa„  sinP^  AP  + V^  sina  cosp„  Aq 

e e e e e e ^ 

- V sinP„  Ar) 

e e 


Yg  AP  + Yg  A6a  + Y5  A <5 r + mg  cosGe  cos<pe  A 4> 
a r 


mg  sin0e  sin4>e  A0 
in 

- Ve  sinae  cosBe  Ap) 


m(sinP^  AV  + V cosp  AP  + V cosa  cosp  Ar 

G G G 0 0 0 


(79) 


(80) 
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(81) 


Zv  AV  + Z ^ Aa  + Zq  Aq  + Z ^ ASg  - mg  cosOe  sin^e  A»J> 

- mg  sinOe  cosife  AO 

■ m(sina  cosp  AV  + V cosaQ  cosp  An 

G G G G G 

- V sina  sinp  AP  + V sinP  Ap 

G G G G g 

- V cosa  cosp  Aq) 

e e e 

La  Aa  + Lg  Ap  + L Ap  + Ar  + AS  + Lg  AS 

. P . a r 

= Ix  Ap  - Izx  Ar  (82) 


M Aa  + M Aq  + M AS 
a 9 5_ 

= Iy  Aq 

N Aa  + N AP  + N Ap  + N Ar  + N AS  + N AS 
a P P _ r Aaafirr 

■ 4r  - !zx 

A<J>  = Ap  + sin>t>e  tanOe  Aq  + cos^e  tanOg  Ar 
A0  = cosi)>e  Aq  - sinie  Ar 


(83) 


(84) 

(85) 

(86) 


Equations  (79)  through  (86)  are  then  non-dimensionalized 
using  Reference  2,  except  for  AP  terms.  As  in  Hawkins' 
study,  to  avoid  confusion  with  the  coupled  equations,  the 
following  relation  is  used: 

AP  - c Ap  (87) 

Thus,  all  terms  are  non-dimensionalized  in  time  by  the  same 
factor.  The  non-dimensional  equations,  after  combining  terms 
and  rearranging  become : 

A A 

- 2y  cosag  cosPg  AV  + 2 u sinae  cosPe  Aa 

/N 

+ 2\i  cos ag  sinPe  Ap 
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- 2CX  AV  - Cx  Aa  + (2m  sinae  cosBe  - Cx  )Aq 

a a q 

+ CT,  cos 0 A0  - 2u  sinO  Ar  - C A 6 

we  e M G X6  s 


(88) 


- 2m  sinBe  AV  - 2m  cosBe  Ap 

= C,,  sinO  sin<t>  AO  - Cv  AP  + 2m  cosa  cosp  Ar 
we  6 e yB  M e 6 

- 2^i  sinae  cosBe  Ap  - Cw  cosOe  cos|e  A4> 

AR  e 


- C A6a  - C A6r 

‘^a  La 


(89) 


- 2m  sinae  cosBe  AV  - 2m  cosae  cos3e  Aa 

A 

+ 2m  sinae  sinPe  AO 

= - 2C  AV  - C_  Aa  - (2m  cosa  cosp  + C )Aq 
Z a q 

+ C sinO  cosi|i  AO  + 2m  sinOe  Ap 


AR 


+ Cw  cosOe  sirM(>e  Aif>  - Cz  AiSg 
e 6 


(90) 


AR  I Ar  - AR  Ix  Ap 

* A A A> 

= - C i Aa  - Ci  AO  - Ar  - Ap 

a a 8 r p 

- Cl  A6r 


C1  A*a 
1 5 

a 


(91) 


- Iy  M - - C ' Cm  Aq  - C A5 

y a q A 


(92) 


- AR  Iz  Ar  + AR  Izx  Ap 

" - Cna  A*  - % A*  - Cnr  A*  - Cnp  AP 

- C AAa  - Cn,  AAr 
6a  6r 


A$  = sinif e tanO  Aq  + co s^e  tanOg  Ar  + 1 Ap 


AR 


AR 


(93) 

(94) 


(95) 


A6  - cos|e  Aq  - sin*f>e  Ar 

AR 

Equations  (88)  through  (95)  are  eight  non-dimensional  linear 

A A A A A 

equations  with  the  eight  unknowns  being  AV,  Am , Aq , AO,  AB, 

A A A 

Ar,  Ap,  and  A*J» . Analysis  of  this  set  of  equations  alone  is 
not  sufficient  since  the  YF-16  is  statically  unstable.  There- 
fore, the  equations  have  to  be  supplemented  with  the  aircraft 
flight  control  system. 

Development  of  the  Longitudinal  Augmentation  Equations 

Appendix  B is  the  diagram  for  the  entire  FCS  for  the 
YF-16.  It  was  desired,  however,  to  include  for  the  stability 
analysis  purposes  only  that  portion  of  the  FCS  which  was 
essential.  The  longitudinal  FCS  was  required  to  give  the 
aircraft  positive  static  stability  and  was,  therefore,  in- 
cluded in  the  linear  analysis.  A simplified  diagram  of  the 
longitudinal  FCS  is  shown  in  Fig.  1.  The  system  reduces  to 
the  following  differential  equation: 

(*)  (s)  (O 

6 + 61  6 + 1385  5 + 14075  A + 57750  6 + 450 00  A 

s s s s s s 

(O  (O  

= b ot  + b a + b ot  + b u + b a+b  a 

1 2 3 4 5 6 

(s)  (4)  

+ c q + c q + c q + c q + c q 

1 2 3 4 5 

(5)  (4) 

+ (1125  F}  Ve/g)  0 + (23625  F?  Ve/g)  0 

+ (163125  Fa  Ve/g)  ’0  + (421875  F Ve/g)  V 

+ (281250  F s Ve/g)  0 (96) 


16 


The  values  of  the  coefficients  for  a and  q vary  according  to 
the  values  of  a as  follows: 

For  a < 15: 


b 

1 

- 1125  F 

Vs 

(97) 

b 

2 

s 

200  F - 

4 

23625  F Ve/g 

3 c 

(98) 

b 

3 

= 

6200  F 

4 

- 163125  F V / g 

3 e 

(99) 

b 

4 

= 

51000  F 

4 

- 421875  F VQ/g 

3 e 

(100) 

b 

5 

= 

45000  F 

4 

- 281250  F V /g 
s e & 

(101) 

b 

6 

= 

0 

(102) 

C 

1 

s 

30  F 

3 

(103) 

c 

2 

=x 

1050  F 

3 

(104) 

C 

3 

= 

12750  F 

3 

(105) 

C 

4 

= 

63750  F 

3 

(106) 

c 

5 

= 

112500  F 

3 

(107) 

O 

15 

< a < 20.4°: 

b 

j 

= 

- 1125  F 

V /g 

3 e 6 

(108) 

b 

2 

22 

241.5  F 

3 

+ 200  F - 23625  F Vj g 

4 3 e 

(109) 

b 

3 

= 

6279  F 

3 

+ 6200  F - 163125  F V/g 

4 3 e 

(110) 

b 

4 

22 

48300  F 

3 

+ 51000  F - 421875  F V^/g 

4 3 e 

(111) 

b 

5 

S3 

132825  F 

+ 45000  F - 281250  F V / g 

3 4 3 e 

(112) 

b 

6 

= 

9056.5  F 

3 

(113) 

c 

1 

<= 

16.9  F2 

3 

+ 30  F 

3 

(114) 

c 

2 

591.5  F2 

3 

+ 1050  F 

3 

(115) 

C 

3 

= 

7182.5  F 

2 + 12750  F 

3 3 

(116) 

C 

4 

= 

35912.5 

F2  + 63750  F 

3 3 

(117) 
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i 

; 

I 

i 

! 

! 

; 

i 

! 

: 

i 


c 

= 

63375  F2  + 112500  F 

(118) 

5 

3 3 

o 

For  a > 

20 

1.4: 

b 

i 

= 

- 1125  F V /g 

3 c 

(119) 

b 

2 

= 

491.5  F + 200  F - 23625  F V / g 

3 4 3 e 

(120) 

b 

3 

= 

15279  F + 6200  F - 163125  F V /g 

3 4 3 e 

(121) 

b 

= 

150800  F + 51000  F - 421875  F V /g 

(122) 

4 

3 4 3 e 

b 

= 

507825  F + 45000  F - 281250  F V /g 

(123) 

5 

3 4 3 e 

b 

6 

= 

371812.5  F 

3 

(124) 

C 

s 

34.4  F2  + 30  F 

(125) 

1 

3 3 

c 

= 

1379  F2  + 1050  F 

(126) 

2 

3 3 

C 

= 

19869.5  F2  + 12750  F 

(127) 

3 

3 3 

C 

= 

121225  F2  + 63750  F 

(128) 

4 

3 3 

C 

s 

= 

260250  F2  + 112500  F 

3 3 

(129) 

F and  F are  variable  gains.  (See  Appendix  B and  Appendix  C.) 

3 4 

Using  Ogata  (Ref  5:675)  this  differential  equation  was 
put  in  state  space  form.  Since  the  forcing  functions  (a,  q, 
and  9)  involve  derivative  terms,  the  form  that  guarantees  a 
unique  solution  is  the  last  six  equations  of  the  set  of  equa- 
tions described  by  equation  (130) . 


Stability  Analysis 

A stability  analysis  of  the  equilibrium  points  obtained 
from  the  trim  program  was  performed  using  the  linear  equations 
coupled  with  the  longitudinal  FCS . Equations  (88)  through 
(95)  were  reordered  and  written  in  the  following  state  space 


form: 


[a][x]  « [b][x]  + [C][U] 
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(130) 


where 


[B]  = 


B ' B 

1 1 1 12 

B i B 

L 2 1 I 2 2. 


(138) 


[»  ]' 
*-  1 1J 


[b  ], 
L 2 1J 


[b  ]. 
L 2 2J 
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b 

b 

b 

b 

b 

0 

0 

1 1 

1 2 

1 3 

1 9 

1 5 

1 6 

b 

b 

b 

b 

0 

b 

b 

2 1 

2 2 

2 3 

2 9 

2 5 

2 7 

2 8 

0 

b 

b 

0 

b 

0 

0 

0 

3 2 

3 3 

3 5 

0 

0 

b 

0 

0 

b 

0 

0 

9 3 

9 6 

b 

0 

b 

b 

b 

b 

b 

S 1 

5 2 

5 9 

5 5 

5 6 

5 7 

5 8 

0 

b 

0 

0 

b 

b 

b 

0 

6 2 

6 5 

6 6 

6 7 

0 

b 

0 

0 

b 

b 

b 

0 

? 2 

7 5 

7 6 

7 7 

0 

0 

b 

0 

0 

b 

b 

0 

9 3 

8 6 

8 7 

• 

0] 

. 

0 

b 

b 

b 

0 

0 

0 

0 

92 

93 

99 

0 

b 

b 

b 

0 

0 

0 

0 

102 

103 

109 

0 

b 

b 

b 

0 

0 

0 

0 

112 

113 

119 

0 

b 

b 

b 

0 

0 

0 

0 

122 

123 

129 

0 

b 

b 

b 

0 

0 

0 

0 

132 

133 

139 

0 

b 

b 

b 

0 

0 

0 

0 

- 

192 

193 

199 

- 

"o 

b 

0 

0 

0 

0 ~ 

910 

0 

0 

b 

0 

0 

0 

1011 

0 

0 

0 

b 

0 

0 

1112 

0 

0 

0 

0 

b 

0 

1213 

0 

0 

0 
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(140) 

(141) 


(142) 


0 


(143) 


[C]« 


0 0 
0 0 
0 0 


c c 

SI  5 i 

c c 

6 1 6 2 

C C 

7 1 7 2 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

0 0 


The  elements  for  matrices  [A]  and  [R]  are  listed  in  Appendix 
D.  Note  that  the  last  six  rows  of  matrix  [r]  vary  for  the 
range  of  a.  The  [c]  matrix  terms  are  valid  for  all  values 
of  a . 


r 


Equation  (130)  is 

[*Ma-'][bKI 


rewritten  as: 

+ [a*  as m 


(144) 

(145) 

(146) 

(147) 

(148) 
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The  stability  of  each  equilibrium  position  was  checked  with 
no  control  inputs.  Therefore,  [if]  = 0,  and  equation  (148) 
becomes : 

[X]  = [A*  0[B]  [X]  (149) 

The  eigenvalue  problem  to  be  solved  is  then: 

[[A*  TO  - X [i]]  [X]=  0 (150) 

which  leads  to  the  condition: 

Det  [[A"®  " X[l]]  =0  (151) 

This  solution  is  a polynomial  of  degree  14: 


14  + 

m 

X1  3 + 

m X 1 2 

+ m X 1 1 + 

m X 1 0 + 

m X 9 

i 

2 

3 

4 

5 

+ m 

X8 

+ m X 

7 + m 

X6  + m X5 

+ m X“ 

+ m X 3 

6 

7 

8 

9 

1 0 

i l 

+ m 

X 

2 + m 

X + m 

= 0 

(152) 

1 2 

1 3 

1 4 

For  the  system  to  be  stable,  all  14  roots  had  to  be  negative. 
Eq.  (152)  was  solved  using  the  computer  program  in  Appendix 
D.  In  addition,  the  eigenvectors  for  the  associated  eigen- 
values were  calculated  for  use  in  identifying  the  eight  roots 
traditionally  used  to  describe  aircraft  characteristics.  The 
six  other  roots  were  due  to  the  flight  control  system  and 
were  of  no  interest. 
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III.  Non-linear  Analysis 


Non-linear  Computer  Program 

The  computer  program  used  for  the  non-linear  analysis 
(Appendix  D)  was  developed  by  the  Aircraft  Group  of  the  Air 
Force  Flight  Dynamics  Labratorv.  It  was  modified  for  this 
study  so  that  the  updated  data  package  could  be  used.  The 
program  contained  the  entire  FCS  as  shown  in  Appendix  R. 

The  inputs  for  the  program  were  designed  to  be  pilot  forces 
for  aileron,  rudder,  and  elevator.  However,  for  purposes 
of  determining  perturbation  effects,  the  program  was  modified 
so  that  perturbations  of  one-half  second  could  be  artificially 
introduced  once  the  aircraft  was  trimmed.  Control  forces 
were  held  constant  once  trimmed  conditions  were  achieved. 

Trim  times  of  15  seconds  and  20  seconds  were  used.  Although 
maintaining  constant  trim  forces  was  somewhat  of  an  artifi- 
cial situation,  this  approach  was  necessary  to  isolate  the 
aircraft's  behavior. 

To  verify  the  accuracy  of  both  the  data  and  the  com- 
puter model,  several  simulations  were  made  with  control 
force  pulses  and  step  inputs.  The  simulations  matched  the 
corresponding  maneuver  from  the  flight  tests  very  closely. 

Perturbed  Motion  Analysis 

During  this  portion  of  the  study,  emphasis  was  on  what 
severe  aircraft  motions  were  encountered  and  if  recovery 
from  the  motions  could  have  been  accomplished.  In  the  inter- 
est of  computer  resources,  it  was  decided  to  use  only  one 
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lateral  perturbation.  Therefore,  yaw  perturbations  were  used 
as  well  as  pitch  perturbations.  Only  one  perturbation  was 
introduced  per  simulation  in  order  to  try  to  separate  what 
motions  could  be  expected  from  longitudinal  perturbations  and 
from  lateral  perturbations.  For  the  case  of  yaw  perturbations, 
a simple  controller  was  modelled  to  control  yaw  divergence. 
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IV.  Results 


Alpha-Beta  Envelope 

The  a-B  envelope  for  which  the  aircraft  could  be  trimmed 
at  sea  level  for  rectilinear  flight  using  the  relationships 
derived  from  eqs . (22)  through  (27)  is  depicted  in  Fig.  2. 
This  envelope  does  not  account  for  restricted  control  deflec- 
tions due  to  the  ARI . It  merely  defines  where  sufficient 
control  authority  is  available. 

Within  the  a- 6 control  envelope  the  linearly-derived 
stability  envelope  is  also  shown  in  Fig.  2.  The  eigenvalue 
computer  program  provided  14  eigenvalues  for  each  given  B 
and  a.  Of  these  values  eight  represented  traditional  air- 
craft behavior.  Except  at  large  sideslip  angles,  the  short 
period,  phugoid,  dutch  roll,  spiral  mode,  and  roll  mode  were 
easily  distinguished.  Values  of  the  dutch  roll  natural  fre- 
quencies calculated  from  the  eigenvalues  agreed  well  with 
flight  test  data.  The  damping  ratios  were  approximately 
one-half  of  flight  test  values  for  ct  up  to  about  20  degrees. 
However,  these  damping  ratios  were  within  Category  B and  C, 
Level  1 handling  qualities.  The  region  below  a = 7° was  un- 
stable  largely  due  to  the  spiral  mode.  Time  to  double  ampli- 
tude (T  ),  however,  was  well  within  Category  A,  Level  1. 

2 

Discrepencies  in  lateral  characteristics  were  expected  since 
the  lateral  FCS  was  not  modelled  in  the  eigenvalue  problem. 
Since  no  lateral  instabilities  of  any  importance  were  dis- 
covered, not  including  the  lateral  FCS  seemed  justified. 

Instability  in  the  phugoid  mode  primarily  accounted  for 
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the  portion  of  the  unstable  region  above  a = 7°.  The  only 
area  where  Level  3 handling  qualities  were  not  satisfied 
was  in  the  area  noted  in  Fig.  2.  In  this  region  T was 
calculated  to  be  as  low  as  ten  seconds.  The  only  region 
that  was  of  concern  for  this  study  was  above  approximately 
a « 10°. 

The  linearized  a-8  control  boundary  and  stability  bound- 
ary provided  guidelines  as  to  where  the  non-linear  computer 
program  should  give  stable  trim  points.  Information  was 
also  gained  as  to  where  in  the  boundary  one  or  more  control 
surfaces  should  become  saturated. 

Perturbed  Motion  Analysis 

Points  in  the  a- 6 boundary  between  a = 10'  and  a - 25 
in  increments  of  two  degrees  were  tested  for  dynamic  sta- 
bility with  various  perturbations.  Fig.  3 shows  regions 
where  the  non-linear  computer  program  was  not  able  to  trim 
the  aircraft.  Analysis  of  the  lower  region  showed  £ was 
maintained,  but  u wandered  to  the  upper  region  and  never 
decreased  in  value.  This  supported  the  linear  analysis  of 
an  unstable  phugoid  region.  The  upper  region  seemed  to  be 
limited  by  control  surface  deflection.  Some  restriction  in 
the  control  boundary  was  anticipated  since  the  ARI  was  part 
of  the  non-linear  model.  Except  for  the  two  regions  noted, 
the  rest  of  the  a-S  control  boundary  above  a = 10  predicted 
accurately  where  trimmed  flight  was  possible. 

As  mentioned  in  Section  III,  only  pitch  and  yaw  per- 
turbations were  introduced.  The  size  of  perturbations  that 
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caused  uncontrolled  aircraft  motion  were  20  deg/sec  and 
30  deg/sec.  These  were  rates  that  could  be  generated  by  the 
pilot  as  he  begins  a variety  of  maneuvers  from  trimmed  flight. 
Lesser  perturbations  resulted  in  at  worst  roll  rates  of 
40  deg/sec  or  merely  a negative  0 with  a and  6 remaining 
constant.  In  all  cases  with  perturbation  rates  below  20 
deg/sec,  no  control  surfaces  were  saturated,  and  recovery 
would  have  been  easily  obtained. 

The  areas  shown  in  Figs.  4 and  5 for  pitch  perturbations 
and  in  Figs.  7 and  8 for  yaw  perturbations  represent  equi- 
librium conditions  from  which  the  aircraft  was  perturbed  and 
the  resulting  motion  was  highly  uncontrollable.  Uncontrol- 
lable for  purposes  of  this  study  means  a condition  where  the 
resulting  motion  was  very  erratic  and  disorienting,  with  one 
or  more  control  surfaces  reaching  the  limits  of  travel  and 
not  affecting  a recovery. 

Pitch  Rate  Perturbations 

The  pitch  rate  perturbations  caused  four  types  of  uncon- 
trolled motion.  The  types  of  motion  and  the  average  behavior 
of  the  aircraft  for  these  motions  are  listed  below: 

(1)  Erect  Spin  (Appendix  E) : 
r = - 60  deg/sec 
0 oscillated  around  - 70  deg 
<j>  within  ±30  deg  of  wings  level 
q highly  oscillatory  and  ranged  between  3.6 


and  4.2  sec/cycle 


a diverged  oscillating  around  60  deg 


V 
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(2)  Inverted  Spin  (Appendix  G) : 
r « -50  to  -70  deg/sec 

0 oscillated  around  -70  deg 
ip  within  +20  deg  of  -150  deg 
q oscillated  at  3.0  to  4.5  sec/cycle 
a oscillated  around  -50  deg 

(3)  Rolling  Departure  (Appendix  I) : 
p = 160  deg/sec 

r oscillated  from  -20  to  -60  deg/sec 
0 oscillated  around  -70  deg 
a,g  diverged  to  negative  values 
q oscillated  at  2.5  sec/cycle 

(4)  Rolling  Departure  to  an  Inverted  Spin  (Appendix  J) : 
r diverged  beyond  -110  deg/sec 

0 oscillated  around  -20  deg 
p within  ±20  deg  of  -160  deg 
q oscillated  at  2.5  sec/cycle 
a oscillated  around  -20  deg  while  rolling  and 
oscillated  around  -50  deg  while  inverted. 


Of  the  four  types  of  departures  due  to  pitch  perturba- 
tions, the  rolling  departures  were  the  most  abrupt  motions. 

The  inverted  spins  that  developed  from  rolling  departures 
had  yaw  rates  three  times  greater  than  yaw  rates  of  spins 
encountered  from  any  other  type  entry.  In  addition,  the  I 


pitch  oscillations  were  twice  those  from  other  spins.  Var- 
iations between  q = 30  deg/sec  and  q = 20  deg/sec  had  little 
effect  on  the  resulting  motion.  Both  pitch  rates  produced 
the  same  four  types  of  motion.  Each  rate  had  a slightly 
different  region  of  influence  on  the  a-6  boundary  (Figs.  4 
and  5).  Fig.  6 shows  the  lines  which  separate  where  the 
controllable  motions  and  the  uncontrollable  motions  occurred 
for  both  pitch  rates.  The  cross-hatched  regions  are  of  inter- 
est because  these  regions  are  where  q = 20  deg/sec  caused 
uncontrolled  motion  while  q = 30  deg/sec  caused  only  small  roll 
rates  (30  to  40  deg/sec).  Appendix  E and  Appendix  K compare 
a trim  point  where  this  occurred. 

Yaw  Rate  Perturbations 

Two  classes  of  departures  resulted  from  yaw  perturbations, 
erect  spins  and  erect  spins  that  transitioned  to  inverted 
spins . 


(1)  Erect  Spin  (Appendix  F) : 
r = -40  to  -80  deg/sec 

0 oscillated  around  -30  deg 
within  ±30  deg  of  wings  level 
q oscillated  at  3.0  to  4.5  sec/cycle 
a diverged  oscillating  around  60  deg 

(2)  Erect  Spin  to  Inverted  Spin  (Appendix  H) 
While  erect,  same  as  above 

While  inverted: 


r diverged  beyond  -50  deg/ sec 
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Fig.  4 Departure  Region  for  a = 20  deg/sec 
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0 oscillated  around  -40  deg 
<p  within  .f20  deg  of  -160  deg 
q oscillated  at  4 sec/cycle 
a oscillated  around -50  deg  while  inverted 

Varying  yaw  perturbations  produced  the  same  sort  of 
results  as  varying  pitch  perturbations.  The  a- 6 region 
influenced  by  r =-30  deg/sec  was  about  twice  that  for 
r =-20  deg/sec.  Fig.  9 shows  a region  where  the  lesser 
rate  caused  uncontrolled  motion  while  the  larger  rate 
resulted  in  only  slow  rolls  (30  to  50  deg/sec).  Appendix 
F and  Appendix  K.  compare  a point  where  this  occurred. 

Comparison  of  the  uncontrolled  motions  produced  by  pitch 
and  yaw  perturbations  revealed  five  basic  types  of  motion: 

(1)  Erect  Spin 

(2)  Erect  Spin  transitioning  to  Inverted  Spin  (yaw 
perturbations  only) 

(3)  Inverted  Spin  (pitch  perturbations  only) 

(4)  Rolling  Departure  (pitch  perturbations  only) 

(5)  Rolling  Departure  to  Inverted  Spin  (pitch  pertur- 
bations only) 

Regardless  of  the  type  perturbation  that  caused  the  out  of 
control  motion,  all  erect  spins  and  all  inverted  spins  were 
very  similar  in  behavior. 

As  a basis  for  comparing  the  erect  spins  described  by 
this  study,  the  one  actual  erect  spin  that  the  YF-16  encoun- 
tered during  its  flight  tests  is  described  below  (Ref  1:234): 
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r = -10  to  -50  deg/sec 
0 = -75  deg 

$ within  ±30  deg  of  wings  level 

q the  primary  oscillatory  motion  at  4.2  sec/cycle 
3 oscillated  from  -25  to  +25  deg 
a was  above  30  deg 

Control  Surfaces  Analysis 

Analysis  of  the  aircraft's  uncontrolled  motions  leads 
to  three  modes  that  could  be  prevented:  erect  spins,  in- 
verted spins,  and  rolling  departures.  Average  control  sur- 
face deflections  during  these  motions  were  as  follows: 

(1)  Erect  Spin  (left  turns) : 

6r  oscillated  0 to  25  deg 

5 =25  deg 
s 

6 oscillated  -5  to  -20  deg 

a ° 

(2)  Inverted  Spin  (right  turns) : 

6^  = 25  to  30  deg 

6 = ±5  deg  of  neutral 

s 

6 = 5 to  -20  deg 

a 

(3)  Rolling  Departure  (right  rolls) : 

6 =30  deg 

r 

6 - -25  deg 

s 

6 = -10  to  10  deg 

a 

The  stabilator  in  all  cases  of  erect  spins  was  fully  de- 
flected in  order  to  try  to  reduce  a.  However,  low  dynamic 
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pressure  apparently  prevented  the  stabilator  from  being 
effective.  This  was  the  case  in  the  actual  departures  and 
the  spin  which  the  aircraft  encountered  in  flight  tests.  For 
the  erect  and  inverted  spins,  rudder  was  being  held  in  the 
direction  of  the  spin.  A controller  that  commanded  at  least 
neutral  rudder,  or  possibly  opposite  rudder  once  the  spin 
has  developed  would  be  seemingly  warranted. 

The  rolling  departures  appear  to  be  the  most  difficult 
departures  to  control.  The  rudder  is  already  deflected  so 
as  to  oppose  the  resulting  roll  direction.  Sufficient  aile- 
rons seem  to  be  available  to  stop  the  roll,  however.  Rolling 
departures  resulted  exclusively  from  pitch  perturbations  and 
occurred  within  two  degrees  of  the  right  side  of  the  a-F 
control  boundary  for  a < 21  deg.  It  appears  that  a controller 
which  prevents  the  aircraft  from  entering  that  region  would 
be  warranted. 

Yaw  Departure  Control  1 e r 

A simple  yaw  departure  controller  was  modelled  which 

used  a,  6,  and  yaw  rate  as  its  inputs.  The  regions  of  Figs. 

7 and  8 were  written  for  8 as  a function  of  a.  When  yaw 

rates  greater  than  20  deg/sec  or  30  deg/sec  were  generated 

and  the  aircraft  was  in  the  appropriate  u-F  region,  rudder 

I 

was  commanded  to  counteract  any  increase  in  8.  The  controller 
was  effective  in  keeping  the  aircraft  from  entering  a spin 
(Appendix  N) . The  FCS  was  able  to  control  a once  the  lat- 
eral motion  dampened  out.  In  all  cases  the  controller  was 
needed  for  less  than  one-half  second. 
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V.  Conclusions  and  Recommendations 

Conclusions 

The  results  of  pitch  and  yaw  perturbations  analysis 
provide  the  following  conclusions  based  on  the  YF-16  com- 
puter model  as  described  in  this  study: 

(1)  The  aircraft  exhibits  a high  degree  of  suscepti- 
bility to  loss  of  control  when  operating  at  high 
angles  of  attack  in  asymmetric  flight  when  exposed 
to  moderate  pitch  and  yaw  perturbations. 

(2)  The  uncontrolled  motions  fall  into  three  basic 

categories:  erect  spins,  inverted  spins,  and 

rolling  departures. 

(3)  Sufficient  control  authority  seems  to  be  available 
to  stop  the  primary  undesirable  motions  of  yaw  and 
roll. 

(4)  A departure  controller  that  monitors  a,  3 and 
changes  in  pitch  rate  and  yaw  rate  appears  to  be 
warranted  if  the  aircraft  is  to  be  allowed  to  ma- 
neuver with  full  control  surface  authority. 

Recommendations 

This  study  should  be  extended  in  the  following  ways: 

(1)  A departure  controller  should  be  investigated  for 
incorporation  into  the  FCS.  Once  the  obvious 
undesirable  aircraft  motions  are  controlled,  an 
analysis  of  any  limitations  imposed  by  the  con- 
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troller  should  be  made. 

(2)  The  aircraft  should  be  analyzed  at  various  alti- 
tudes to  determine  if  the  sea  level  a-g  boundary 
and  corresponding  areas  of  uncontrolled  motion  are 
representative  of  all  altitudes. 

(3)  Asymmetrical  flight  should  be  investigated  in  the 
landing  configuration.  This  area  provides  the 
greatest  safety  threat  since  any  uncontrolled 
motion  would  have  extremely  limited  time  and  alti- 
tude for  recovery. 

(4)  Turning  flight  should  be  investigated  for  diver- 
gence problems. 
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Appendix  A 


YF-16  Physical 

W = 20,860  lb 
S = 283  ft2 
c = 10.94  ft 
b = 29  ft 
AR  = 3.0 


Characteristics 

Ix  = 6,100  slugs-ft2 
Iy.  = 51,400  slugs-ft 
Iz  = 57,000  slugs-ft 
I = 160  slugs-ft2 


Control  Surfaces  Deflection  Limits 
(Primary  Control  Surfaces) 


Aileron  Deflection 

20 

deg 

up 

20 

deg 

down 

Rudder  Deflection 

30 

deg 

left 

30 

deg 

right 

Stabilator  Deflection 

25 

deg 

up 

25 

deg 

down 

Note:  TEL  = Trailing  Edge  Left 
TED  - Trailing  Edge  Down 

6f  “ 6a 
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CALL  SCALF(Q,?.»L00P,1) 

CALL  SCALED, 2., LOOP, 1) 

CALL  STALE <THETA, 2. .LOOP, 1) 
CALL  SrALE(PSI,2»,L00p»l) 
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Appendix  E 
Erect  Spin 

a = 21  deg,  g = 10  deg 
I = 20  deg/sec  at  t = 15  sec 
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Fig.  E-3  0,  ^ , 1(1  vs . time,  a = 21  deg,  £ *>  10  deg, 

q - 20  deg/sec  at  t = 15  sec 
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Appendix  F 
Erect  Spin 

ct  = 21  deg,  g = 10  deg 
r = -20  deg/sec  at  t = 15  sec 
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Fig.  F-2  a,  6.  and  Velocity  vs.  time,  a = 21  deg, 

6 * 10  deg,  r = -20  deg/sec  at  t = 15  sec 
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Fig.  F-4  Control  Forces  and  Deflections  vs.  time,  a = 21  deg, 
0 = 10  deg,  r = -20  deg/sec  at  t = 15  sec 
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Inverted  Spin 
a = 23  deg,  g = 13  deg 
q = 30  deg/sec  at  t = 20  sec 
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Fig.  H-l  p,  q,  and  r vs.  time,  a = 19  deg,  3=4  deg, 
r = -30  deg/sec  at  t = 15  sec 
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Fig.  H-5  N2  and  Altitude  vs.  time,  a = 19  deg,  (5-4  deg, 
r * -30  deg/sec  at  t = 15  sec 


Fig.  1-1  p,  q,  and  r vs.  time,  a = 19  deg,  3=9  deg, 
q = 20  deg/sec  at  t = 15  sec 
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Fig.  1-2  a,  8,  and  Velocity  vs.  time,  a = 19  deg, 
8 ■ 9 deg,  q =*  20  deg/sec  at  t = 15  sec 
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AILERON  FORCE-LBS  ELEVATOR  FORCE-LBS  RUDDER  FORCE-LBS 
-8.00  0.00  8.00  0.00  1.00  2.00  0.00  1.00  2.00 
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a - 17  deg,  B = 9 deg 
= 20  deg/sec  at  t = 15  sec 
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Fig.  J -2  a,  $ , and  Velocity  vs.  time,  ct  = 17  dog, 
0 «■  9 deg,  q = 20  deg/sec  at  t = 15  sec 
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Fig.  K-3  0,  4>,  and  ip  vs.  time,  at  = 21  deg,  & = 10  deg, 

q = 20  deg/sec  at  t = 15  sec 
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Fig.  K-4  Control  Forces  and  Deflections  vs.  time,  a 
6 = 10  deg,  q = 20  seg /sec  at  t = 15  sec 


= 21  deg, 
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Fig.  L-2  a,  0,  and  Velocity  vs.  time,  a = 25  deg, 
6=6  deg,  r = -20  deg/sec  at  t = 15  sec 
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ILERON  FORCE-LBS  ELEVRJOR  FORCE-LBS  RUDDER  FORCE-LBS 
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Fig.  L-4  Control  Forces  and  Del lections  vs.  time,  a 
0-6  deg,  r = -20  deg/sec  at  t - 15  sec 


= deg, 


135 


6 deg. 


Fig.  L-5  Nz  and  Altitude  vs.  time,  .'i  = 25  deg,  £ = 
r = -20  deg/sec  at  t = 15  sec 


Appendix  M 


No  Uncontro 1 led  Motion 
a ■ 25  deg , |<  ° 6 do^ 
r ■"  -30  dep,/sec  nt  L ■=  15  sec 
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Fig.  M-2  a,  0,  and  Velocity  vs.  time,  a = 25  deg, 
8-6  deg,  r - -30  deg/sec  at  t = 15  sec 
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Fig.  M-4  Cor 


— 1 — 
20. 

00 

1 

30.00 

1 — 

40.00 

1 

50.00 

'IM 

E - 

SEC 

Control  Forces  and  Deflections  vs.  time,  a 
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Fig.  N-2  a.  0,  and  Velocity  vs.  time,  ot  = 19  deg, 
P « 4 deg,  r = -30  deg/ sec  at  t = 15  sec 
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His  last  assignment  prior  to  reporting  to  AFIT  in  September 
1976  was  to  Craig  AFB,  Alabama.  While  at  Craig  AFB,  he  was 
a T-38  Instructor  Pilot,  and  his  duties  included  Squadron 
Executive  Officer  for  the  52nd  Flying  Training  Squadron  and 
Wing  Flying  Safety  Officer. 
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